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Abstract 


Planetary boundary layer (PBL) processes, including those influencing the PBL depth, 
control many aspects of weather and climate and accurate models of these processes are 
important for forecasting changes in the future. However, evaluation of model estimates of PBL 
depth are difficult because no consensus on PBL depth definition currently exists and various 
methods for estimating this parameter can give results that differ by hundreds of meters or more. 
In order to facilitate comparisons between the Goddard Earth Observation System (GEOS-5) and 
other modeling and observational systems, seven PBL depth estimation methods are used to 
produce PBL depth climatologies and are evaluated and compared here. All seven methods 
evaluate the same atmosphere so all differences are related solely to the definition chosen. These 
methods depend on the scalar diffusivity, bulk and local Richardson numbers, and the diagnosed 
horizontal turbulent kinetic energy (TKE). Results are aggregated by climate class in order to 
allow broad generalizations. 

The various PBL depth estimations give similar midday results with some exceptions. 
One method based on horizontal turbulent kinetic energy produces deeper PBL depths in the 
winter associated with winter stonns. In warm, moist conditions, the method based on a bulk 
Richardson number gives results that are shallower than those given by the methods based on the 
scalar diffusivity. The impact of turbulence driven by radiative cooling at cloud top is most 
significant during the evening transition and along several regions across the oceans and methods 
sensitive to this cooling produce deeper PBL depths where it is most active. Additionally, 
Richardson number-based methods collapse better at night than methods that depend on the 
scalar diffusivity. This feature potentially affects tracer transport. 
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Introduction 


The planetary boundary layer (PBL) is crucial to surface-atmosphere exchanges of 
momentum, energy, moisture, aerosols, carbon, and other chemical tracers. Accurately modeling 
the PBL character and height in global atmospheric general circulation models (AGCMs) has 
significant implications for climate and weather predictions, and is made difficult due to the lack 
of adequate observations. The depth of this layer, in particular, is important for air quality 
studies since pollutants, aerosols, and carbon dioxide that are emitted near the surface are 
turbulently mixed throughout the layer and the depth reflects the amount and profile of the 
turbulence. Inversion studies, which seek to estimate surface fluxes given surface 
concentrations, are explicitly sensitive to PBL depth errors in atmospheric models [ Gurney et al., 
2002 ]. 

The PBL parameterization is flawed in many numerical weather models and this leads to 
prediction errors [ Beljaars , 1995; J off re et al., 2001]. These flaws, combined with only limited 
measurements, contribute to the difficulty of producing a PBL depth comprehensive climatology, 
despite its importance [ Joffre et al, 2001]. In addition, multiple PBL depth definitions exist and 
these can give different results [ Seidel et al., 2010], complicating PBL depth estimation and 
comparisons between models and observations. 

Several studies have attempted to understand the uncertainty associated with the use of 
different PBL depth definitions and found the result to depend substantially on the method 
chosen. Vogelezang and Holtslag [1996] examined the PBL depth by defining it using both bulk 
and gradient Richardson numbers and found that the choice of Richardson number, the critical 
number chosen, and the inclusion of surface friction impacted the results. White et al. [1999] 
used summertime data near Nashville, Tennessee to compare the PBL depth estimated by wind 
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profilers and an airborne differential absorption lidar (DIAL). They found that the different 
datasets produced similar PBL depths under clear conditions (correlation coefficient of 0.94 and 
mean offset of 37 m), but the agreement degraded under cloudy conditions (correlation 
coefficient of 0.87) with a mean difference of about 150 meters. 

Seibert et al. [2000] found that all PBL height definition schemes had deficiencies under 
certain conditions. They examined 7 different datasets and 10 equations commonly used to 
detennine the PBL depth and concluded that the definitions have to be seen in the context of the 
data used. Likewise, Seidel et al. [2010] tested 7 different PBL depth definition methods on 
radiosonde profiles. Using a single dataset, the estimated PBL depth was found to generally 
differ by several hundred meters, depending on the method used. The use of different methods 
in their study produced differences of more than 1 km and even produced different seasonal 
variations. These differences were statistically significant in practically all comparisons between 
methods, introducing a structural uncertainty of 10-100% of climatological means. Seidel et al. 
[2010] concluded that it is necessary to compare different PBL depth estimates using the same 
method. They suggested the development of multiple climatologies using different definitions 
and the use of the appropriate one for the application desired. 

In the present study, seven different methods to compute the PBL depth, based on vertical 
profiles of the scalar diffusivity, the bulk and local Richardson (Ri) numbers, and the horizontal 
component of the turbulent kinetic energy (TKE), are incorporated into the Goddard Earth 
Observation System (GEOS-5) AGCM [Rienecker et al., 2008; Molod et al., 2012] and 
compared using a single climate simulation. In order to provide insight into implications on the 
regional and global climate scale, results are aggregated using Koppen-Geiger climate classes 
[ Peel et al, 2007]. 
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The purpose of this study is two-fold. First, it analyzes differences among the PBL depth 
definitions within the GEOS-5 model using the same atmosphere. Second, it enables a 
comparison among GEOS-5, observations, and other numerical simulations using the same 
definition so as to evaluate the atmospheric profiles rather than the method used. 

The following section provides a model description and a description of the PBL depth 
diagnostics used. The third section presents results, and the final section contains the 
conclusions. 

2. Model and PBL Diagnostics 

2. 1 GEOS-5 Model Description 

The GEOS-5 AGCM is part of the GEOS-5 data assimilation system, an earlier version 
of which was used for the Modem-Era Retrospective Analysis for Research and Applications 
(MERRA) [Rienecker et al., 2011]. The latitude-longitude hydrodynamical core of GEOS-5 
uses the finite volume dynamical core of Lin [2004] and the cubed sphere version is based on 
Putman and Lin [2007]. GEOS-5 includes moist physics with prognostic clouds [Bacmeister et 
al., 2006]. The convective scheme is a modified version of the Relaxed Arakawa-Schubert of 
Moorthi and Suarez [1992], the shortwave radiation scheme is that of Chou and Suarez [1999], 
and Chou et al. [2001] describe the longwave radiation scheme. The Catchment Land Surface 
Model is used to determine fluxes at the land/atmosphere interface [ Poster et al., 2000] and the 
surface layer is determined as in Helfand and Schubert [1995]. The model uses 72 vertical 
pressure layers that transition from terrain following near the surface to pure pressure levels 
above 180 hPa [Rienecker et al, 2008; Molod et al, 2012] and this study uses approximately Vi 
degree horizontal resolution on the cubed sphere. The simulation covers January 1990 through 
May 2013. 


5 


116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 


GEOS-5 includes two atmospheric boundary layer turbulent mixing schemes [Rienecker 
et al., 2008]. The scheme of Louis et al. [1982] is used in conjunction with the scheme of Lock 
et al. [2000]. The Lock scheme is a non-local first order scheme in which the diffusivities are 
computed based on the buoyancy associated with the surface based (positive buoyancy) and the 
cloud-based radiative cooling (negative buoyancy) “plumes”. This scheme has been extended in 
GEOS-5 so that the unstable surface plume calculation includes moist heating and entrainment 
[Rienecker et al, 2008]. The Louis scheme is a first order local scheme in which the turbulence 
diffusivities are computed as functions of the gradient Richardson number. The turbulent length 
scale is assumed to be related to the PBL height as diagnosed from the Louis-Lock combined 
eddy diffusivities. The eddy diffusivities used for the AGCM turbulent diffusion are the larger 
of the Lock or Louis diffusivities at any time step [ Molod et al., 2012]. 

2.2 PBL Depth Diagnostics 

Seven different methods for determining the PBL depth are evaluated using the GEOS-5 
model based on several different output variables (Table 1). All methods evaluate the same 
atmospheric profiles and all differences are related solely to the choice of definition. The PBL 
depth based on the total scalar diffusivity using method 1, is used to compute the turbulent length 
scale [B/ackadar, 1962] for the Louis scheme. This PBL depth is tenned “active” since it feeds 
back to the turbulence scheme and detennines the vertical extent of mixing. This first method 
estimates the PBL depth as the model level above which the scalar diffusivity falls below a 
threshold value of 2 m s' . 

The use of the scalar diffusivity to define the PBL depth is further investigated by using a 
threshold of 10% of the column maximum and linearly interpolating between levels to determine 
the PBL depth. Method 2 uses the total scalar diffusivity and method 3 uses the surface 
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buoyancy driven scalar diffusivity, but neglects the radiatively driven component defined by the 
Lock et al. [2000] scheme. 


The PBL depth definition used by Seidel et al. [2012] and based on the work by 

Vogelezang and Holtslag [1996] is used as the fourth method. Seidel et al. [2012] selected this 

method because of its applicability to radiosondes and model simulations and its suitability for 

convective and stable boundary layers. This method uses a bulk Richardson number given by: 

(^-)(»*» )(*-**) 

U VZ 


fli'bO) = 




Rib is the Richardson number, g is gravity, 0 V is the virtual potential temperature, u and v are the 
horizontal wind components, and z is height. The subscript s denotes the surface and the surface 
winds are assumed zero. The Richardson number is evaluated between the surface and 
successively higher heights, identifying the PBL top as the level at which Ri exceeds a critical 
value of 0.25. The PBL height is found by linearly interpolating between model levels. 
Additionally, two other methods use a different version of the Richardson number evaluated 
between two consecutive levels rather than between the surface and the current height. These 
methods use a local Richardson number calculated as: 

Ri Z (u rl -u„) a +(v xl -v„) a 


Here, z2 represents the level height below the height zl and 0 V without a subscript is the average 
virtual potential temperature between heights zl and z2. The fifth method tests a critical 
Richardson number value of 0 and the sixth tests a critical value of 0.2. 

Finally, since the PBL is generally considered turbulent with only sporadic clear air 
turbulence aloft \Stull, 1988], the TKE due to shear is used to estimate the PBL depth in the 
seventh method. The horizontal TKE is estimated from the wind shear and momentum 
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diffusivity and then the top of the PBL is taken to be the height at which the value falls below a 
threshold value of 10% of the column maximum. 

2.3 Climate Classes 

Peel et al. [2007] recently updated the Koppen-Geiger climate classification, taking 
advantage of advances in data availability and computing power. They did this by using monthly 
mean precipitation and temperature data from over 4000 stations (plus additional data from 
stations reporting only temperature or only precipitation) and interpolating between them using a 
two-dimensional thin-plate spline with tension. The final map is generated on a 0.1°x0.1° grid. 
The highest station density is in the USA, southern Canada, northeast Brazil, Europe, India, 
Japan, and eastern Australia while the lowest station data is located in desert, polar, and some 
tropical regions. 

Peel et al. [2007] used the same classes as the original classification system, but with an 
updated boundary condition between the temperate and cold climate classes. The classification 
consists of five main classes consisting of tropical (A), arid (B), temperate (C), cold (D), and 
polar (E) with further divisions based on seasonal variations in temperature and precipitation. 
Peel et al. [2007] provides a full description of the climate classifications including details on 
how the classification was determined. The broad climate types are relatively insensitive to 
temperature trends, including those from global climate change [ Triantafyllou and Tsonis, 1994; 
Peel et al, 2007] and are intended to represent long tenn mean climate conditions and not year- 
to-year variability. 

In an effort to generalize the results of this analysis, the computed PBL depths are 
aggregated by season onto the Koppen-Geiger climate classes, which are regridded to the 1/2° 
horizontal grid used by GEOS-5 (Figure 1). The Koppen-Geiger climate classes have been used 
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to group rivers worldwide for comparisons of runoff characteristics [ McMahon et al., 1992; Peel 
et al., 2004]. Similarly, Molod and Salmun [2002] successfully used this technique in their study 
investigating the implications of using different land surface modeling approaches. Their study 
aggregated results such as canopy temperature, soil moisture, and turbulent fluxes and they were 
able to use these results to make generalizations that extend to broad climate regions relevant for 
global models. These climate classes are a way to characterize similar remote regions and apply 
findings globally. However, it is important to keep in mind that this classification does not take 
into account other aspects of the climate system relevant to boundary layer processes. 
Differences such as intensity of precipitation, elevation, and overlying subsidence are not 
considered. 

3. Results 

3 . 1 Regimes within Climate Classes 

The Koppen-Geiger climate classes organize remote land regions together based on 
temperature and precipitation. However, these regions can have very different large-scale 
dynamics that also influence the PBL depth such as subsidence, intensity and frequency of 
precipitation, and terrain. Two examples of this type of variability within climate classes are 
presented here. The first example is illustrated in Figure 2, which shows the climatological 
relationship between PBL depth and sensible heat flux for climate class BWh (arid, hot desert) in 
winter. Each point on the scatter plot represents the seasonal mean midday PBL depth and 
sensible heat flux for a GEOS-5 grid cell within the BWh climate class. 

Several different regimes are present associated with different climatic conditions. The 
Australian deserts (the non-black points, colored according to evaporative fraction), in particular, 
have two different regimes associated with them during the winter. The first regime at low 
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2 .... . 

sensible heat flux (<100 W m' ) is associated with evaporative fractions above 0.3 while the 
other regime, characterized by sensible heat fluxes between 100 and 260 W m' and shallower 
PBL depths, has evaporative fractions around 0.2. This second regime has less variation of PBL 
depths with sensible heat, and has a median depth about 260 meters shallower than the first 
regime. This shows that similar climate and physical proximity cannot explain all PBL depth 
variability. 

The second example is illustrated in Figure 3, which shows the relationship between PBL 
depth and 10-meter temperature for the tropical rainforest climate class (Af) colored according to 
relative humidity. In this climate class and the other tropical climate classes, there is a shift in 
the relationship occurring around 302 K. This temperature is near the wilting point for broadleaf 
evergreen, the dominant vegetation type in the tropics. At temperatures above the wilting point, 
the vegetation experiences stress thus severely limiting transpiration and near surface humidity. 
In these drier conditions, less energy goes in to evaporating water and, by energy balance, more 
goes in to sensible heat flux. Since sensible heat is much more efficient at growing the PBL than 
latent heat, the PBL depth increases rapidly with temperature in this drier regime [Avissar and 
Pielke, 1989]. In the regime below the wilting point, transpiration increases with temperature 
and proceeds with little resistance, wetting the lower atmosphere. In this wetter regime, PBL 
depth decreases with temperature. 

These different regimes and sensitivities of PBL depth to different variables must be kept 
in mind when examining climatological boundary layer depth. Although the Koppen-Geiger 
climate classes are able to capture a lot of the variability and are useful for organizing land 
regions in order to make generalizations and simplify the analysis, they do not capture all the 
large-scale climate conditions relevant to boundary layer processes. There will therefore be 
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geographical differences between regions within a climate class that will not be captured by this 
analysis. 

3.2 General Method Behavior 

Generally, the PBL depth definitions produce similar results. Differences are highlighted 
below. Figure 4 shows the seasonal mean diurnal cycle for the cold climate class with dry, hot 
summer (Dsa during summer 4a and winter 4c) and hot, arid desert (BWh during summer 4b and 
winter 4d). For these climate classes, the scalar diffusivity methods using a 10% threshold 
(methods 2 and 3) are insensitive to the use of the radiative plume from the Lock turbulence 
scheme. The PBL depths estimated using the bulk Richardson number (method 4), the TKE 
(method 7), and the three scalar diffusivity methods (methods 1, 2, and 3) give comparable 
midday results over land, giving us confidence in the estimated depth under these conditions. 

Although the horizontal TKE definition gives similar midday results to the scalar 
diffusivity and bulk Richardson number methods under most conditions, during the winter in 
temperate and cold climates, the horizontal TKE method often gives PBL depths that are 100 
meters or more deeper than the other methods (Figure 4c and 4d). This is due to the winter storm 
tracks and associated increase in wind shear. The horizontal TKE is diagnosed in the model by 
multiplying the momentum diffusivity by the local wind shear making the horizontal TKE 
method more sensitive to the wind profile and seasonal changes to it than the other methods. 
The stronger and deeper wind shear thus produces a deeper turbulent layer and therefore a deeper 
diagnosed PBL depth based on this variable. 

The results from the horizontal TKE method, while consistent with other PBL methods 
over land, are inconsistent with the other methods over the oceans. The horizontal TKE derived 
PBL depth gives results as much as 400 meters deeper than the other definitions and have diurnal 
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variability considerably larger. It is therefore suggested that this method not be used to estimate 
marine boundary layer depths. 

The methods based on the local Richardson number estimate PBL depths that are several 
hundred meters shallower at midday than PBL depths using the other methods. This method 
does not depend greatly on the critical Richardson number chosen as the differences between 
PBL depths estimated using a critical value of zero are only slightly shallower than PBL depths 
estimated using a critical value of 0.2. Seidel et al. [2012] similarly found small uncertainty 
associated with the choice of critical value in their study. 

The Richardson number based methods (local and bulk) collapse better at night than the 
definitions based on the scalar diffusivity or TKE. This has implications for estimating the 
shallow nocturnal boundary layer and studies involving tracer transport. For instance, over 
climate class BWh (arid, hot desert representing such areas as the Sahara and Australian deserts, 
Figure 4b), the bulk Richardson number nocturnal PBL is well under 500 meters and the local 
Richardson number PBL depths are only a few hundred meters while the scalar diffusivity 
methods estimate a PBL depth between 1000 and 1500 meters at night during the summer. 

3.3 Bulk Richardson vs. Scalar Diffusivity Methods 

While the bulk Richardson number and scalar diffusivity methods generally give similar 
midday results, under warm, wet conditions, the estimated daily maximum PBL depth found 
using the bulk Richardson number method tends to be shallower than the scalar diffusivity 
methods. This occurs for the tropical rainforest climate class (Figure 5) as well as the other 
tropical climate classes during their rainy seasons and even for temperate climate classes when it 
is both warm and there is a lot of precipitation. This difference means that the Richardson 
number exceeds its critical value at a level below which the scalar diffusivity decreases to its 
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threshold value. This implies a virtual potential temperature inversion occurring within the 
layer of relatively high scalar diffusivity. Figure 6 shows the June, July, and August (JJA) 
seasonal mean vertical profiles of total scalar diffusivity, the scalar diffusivity from the Louis 
parameterization, and the Richardson number from a typical location within the Amazonian 
rainforest. The horizontal dashed lines indicate the PBL depth found using the total scalar 
diffusivity (method 1, Figure 6a) and Richardson number (method 4, Figure 6b). These profiles 
show that the Richardson number becomes stable below where the scalar diffusivity declines. 

This could occur under several different meteorological conditions. There could be a 
turbulent layer aloft that is not fully decoupled from the surface layer that is being detected by 
the scalar diffusivity methods, but not the bulk Richardson number method. If the scalar 
diffusivity predicted by the Louis scheme (Figure 6a) is shallow with its maximum low to the 
ground, it can be expected that the PBL depth found using the bulk Richardson number might 
also be shallow. The Louis turbulence parameterization is dependent upon a local Richardson 
number and so contains some infonnation about the vertical profile of temperature and shear. 
While this is a different form of the Richardson number from what is used in the bulk 
Richardson number method, the Louis scheme can provide information about what to expect 
from the bulk Richardson number method. If, however, the Lock scheme is strongly active aloft 
due to entrainment or radiation, the scalar diffusivity methods may detect a deeper PBL. 

3.4 Impact of radiative plume 

The Lock turbulence scheme extends the nonlocal eddy-viscosity based scheme based on 
Holtslag and Boville [1993] to include the effects of radiative cooling at cloud top by positing a 
turbulent “radiative plume” descending from cloud base due to negative buoyancy associated 
with radiative cooling. In order to test the impact of the associated radiative plume, the scalar 


13 


299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 


diffusivity method using a threshold of 10% of the column maximum was tested both with 
(method 2) and without (method 3) the scalar diffusivity associated with the radiative plume. 
Figure 7 shows the diurnal cycle of the PBL depth difference between the two methods for JJA. 
At all locations, the PBL depth estimated using the radiative plume was deeper or the same 
height as that without the plume. The largest differences occurred over land in the summer 
hemisphere and in the tropics during the evening transition. This result holds for December, 
January, and February (DJF) as well (not shown). The evening timing is because the radiative 
plume is sensitive to cloud top. At night, the scalar diffusivities decrease due to the lack of 
incoming solar radiation, but the diffusivity associated with the radiative plume decreases less 
proportionally to the other diffusivities since the cloud does not dissipate during the evening 
transition. The radiative plume scalar diffusivity thus becomes proportionally more important at 
night and the PBL depth remains deeper. The non-radiative method therefore collapses better at 
night. 

In addition to the maxima over land, there are persistent regions of relatively large 
radiative plume impact over the oceans around 30°N and 45°S. This is due to the way GEOS-5 
estimates the effective liquid cloud particle radius. GEOS-5 uses an empirical model based on 
temperature, pressure, and wind speed to estimate this radius. The large differences over oceans 
are located in regions with a temperature transition, creating a collocation of relatively small 
prescribed effective radii and boundary layer clouds. In cloudy regions, more shortwave 
radiation is absorbed in areas with many small droplets than in areas with fewer large droplets. 
This produces an increase in temperature leading to an enhanced longwave flux, which drives the 
radiative plume. Since the radiative plume is more active in these locations, PBL methods 
sensitive to its impact are deeper than if the radiative plume was neglected. 
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4 . 


Conclusions 


Although the PBL depth is important for AGCMs and has implications for climate and 
weather prediction, observations are limited and no consensus on definition exists. Complicating 
things further, under certain conditions, different definitions can give significantly different 
results. This study examines this issue by using seven different methods to define the PBL depth 
using the same atmosphere so that all differences can be attributed directly to the definition. It 
is, however, important to point out that this is not an extensive validation. Results are 
aggregated to Koppen-Geiger climate classes in order to make broad generalizations and 
simplify the analysis on a global scale. Within these climate classes, different regimes are 
present that are related to differences in temperature and evaporative fraction. These differences 
are not captured when the results are aggregated, but the general behavior is represented. 

Under most conditions, the bulk Richardson number, scalar diffusivity, and horizontal 
TKE methods give similar midday results over land. The horizontal TKE definition is more 
sensitive to winter storms and their associated shear and so estimate deeper midday PBL depths 
during the winter season. The local Richardson number methods are relatively insensitive to the 
critical Richardson number used and estimate PBL depths several hundred meters shallower than 
the other methods. Both the local and bulk Richardson number methods collapse better at night 
than the scalar diffusivity and TKE methods. For instance, over hot, arid deserts, the bulk 
Richardson number method estimates a nocturnal PBL depth that is up to a kilometer shallower 
than the default scalar diffusivity method using a threshold of 2 m s' . 

Under warm, moist conditions, the bulk Richardson number method estimates PBL 
depths that are shallower than those estimated by the scalar diffusivity methods. This indicates 
that the Richardson number is exceeding its threshold value below the level at which the scalar 
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diffusivity decreases to its threshold value. This response is associated with a maximum in the 
Louis scheme scalar diffusivity near the surface and Lock scheme scalar diffusivities with 
maxima higher in the atmosphere. 

The Lock scheme includes the effects of a radiative plume associated with longwave 
cloud base cooling. The impact of this effect on PBL depth was found to have its strongest 
effect over land during the evening transition because of the persistence of cloud cover through 
the diurnal cycle. Additionally, regions of influence were found in the marine boundary layer 
related to the way the effective liquid cloud radius is defined in GEOS-5. 

This study analyzes the sensitivity of the PBL depth in GEOS-5 to the definition chosen 
under different meteorological conditions. Although each definition was evaluated using 
identical atmospheres, under various environments the resulting PBL depth could differ by a 
kilometer or more. This reinforces the recommendation by Seidel et al. [2010] that PBL depths 
be compared only using the same definition. This work also provides multiple climatologies that 
can be used for evaluation and comparisons in the future with other models and observations 
using similar definitions. Since the PBL depth in GEOS-5 is not purely diagnostic, but feeds 
back on the turbulence through the turbulent length scale, the choice of definition is important 
for vertical mixing and tracer transport. This aspect of the model will be evaluated in future 
work. 
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Figure 1 Koppen-Geiger climate classes as determined by Peel et al. [2007] regridded to the 
0.5°x0.5° grid used by GEOS-5. The first letter indicates the broad climate class as tropical (A), 
arid (B), temperate (C), cold (D), and polar (E). 

Figure 2 Scatter plot of PBL depth versus sensible heat flux for the arid, hot desert climate class 
in winter. Each dot represents the seasonal mean midday PBL depth and sensible heat flux. The 
PBL depth is defined using the default scalar diffusivity definition in GEOS-5. The Australian 
deserts are highlighted in color according to evaporative fraction. 

Figure 3 Scatter plot of PBL depth versus 10-meter temperature for the tropical rainforest climate 
class in the winter. Each dot represents the seasonal mean midday PBL depth and sensible heat 
flux. The PBL depth is defined using the default scalar diffusivity definition in GEOS-5. The 
colors highlight the relative humidity. 

Figure 4 Seasonal mean diurnal cycle of PBL depth for climate classes Dsa (Cold with dry, hot 
summers, during summer and winter, 4a and 4c) and BWh (hot, arid desert, during summer and 
winter, 4b and 4d) using 7 different methods for estimating the PBL depth. 

Figure 5 Seasonal mean diurnal cycle of PBL depth for climate classes Af (tropical rainforest) 
during summer (5a) and winter (5b) using 7 different methods for estimating the PBL depth. 
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480 Figure 6 Seasonal mean vertical profile of total and Louis scalar diffusivities (6a) and bulk 

481 Richardson number (6b) for JJA in the Amazonian rainforest. The dashed lines represent the 

482 PBL depth as determined by method 1 (6a) and method 4 (6b). 

483 

484 Figure 7 Diurnal cycle of PBL depth response to radiative plumes during JJA. The figure shows 

485 the scalar diffusivity method using a 10% of the column maximum threshold including the 

486 radiative plume minus the same method, but without the radiative plume. Each subplot is 

487 labeled with the current time in UTC. 
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Figure 1 Koppen-Geiger climate classes as detennined by Peel et al. [2007] regridded to the 
0.5°x0.5° grid used by GEOS-5. The first letter indicates the broad climate class as tropical (A), 
arid (B), temperate (C), cold (D), and polar (E). 
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Hot, Arid Desert in Winter 
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Figure 2 Scatter plot of PBL depth versus sensible heat flux for the arid, hot desert climate class 
in winter. Each dot represents the seasonal mean midday PBL depth and sensible heat flux. The 
PBL depth is defined using the default scalar diffusivity definition in GEOS-5. The Australian 
deserts are highlighted in color according to evaporative fraction. 
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Figure 3 Scatter plot of PBL depth versus 10-meter temperature for the tropical rainforest climate 
class in the winter. Each dot represents the seasonal mean midday PBL depth and sensible heat 
flux. The PBL depth is defined using the default scalar diffusivity definition in GEOS-5. The 
colors highlight the relative humidity. 
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507 Figure 4 Seasonal mean diurnal cycle of PBL depth for climate classes Dsa (Cold with dry, hot 


508 summers, during summer and winter, 4a and 4c) and BWh (hot, arid desert, during summer and 


509 winter, 4b and 4d) using 7 different methods for estimating the PBL depth. 
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Figure 5 Seasonal mean diurnal cycle of PBL depth for climate classes Af (tropical rainforest) 
during summer (5 a) and winter (5b) using 7 different methods for estimating the PBL depth. 
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Figure 6 Seasonal mean vertical profile of total and Louis scalar diffusivities (6a) and bulk 
Richardson number (6b) for JJA in the Amazonian rainforest. The dashed lines represent the 
PBL depth as determined by method 1 (6a) and method 4 (6b). 
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521 Figure 7 Diurnal cycle of PBL depth response to radiative plumes during JJA. The figure shows 

522 the scalar diffusivity method using a 10% of the column maximum threshold including the 

523 radiative plume minus the same method, but without the radiative plume. Each subplot is 

524 labeled with the current time in UTC. 
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525 


526 

527 Summary of PBL depth Methods 


Method 

Abbreviation 

Description 

1 

Kh: 2 threshold 

Uses total scalar diffusivity and a 
threshold of 2 m s' , this is the default 
PBL depth in GEOS-5 

2 

Kh: 10% threshold, rad 

Uses total scalar diffusivity and a 
threshold equal to 10% of the column 
maximum, includes the radiative plume 

3 

Kh: 10% threshold, no rad 

Uses total scalar diffusivity and a 
threshold equal to 10% of the column 
maximum, does not include the radiative 
plume 

4 

Bulk Ri 

Uses the bulk Richardson number used 
by Seidel et al. [2012] and a critical 
value of 0.25 

5 

Ricrit = 0.2 

Uses a Richardson number and a critical 
value of 0.2 

6 

Ricrit = 0 

Uses a Richardson number and a critical 
value of 0 

7 

Horizontal TKE 

Uses the diagnosed horizontal turbulent 
kinetic energy and a threshold of 10% of 
the column maximum 


528 

529 Table 1 List of the PBL depth methods used along with the abbreviations. 
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